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Abstract 

We investigate one-loop weak corrections to the production cross section of two 6-jets at Teva- 
tron and Large Hadron Collider (LHC). We establish that they are small at inclusive level but 
dominant in exclusive observables that have a non-trivial dependence on the helicity structure 
of the hard subprocesses. Such effects can serve as a test of the Standard Model (SM) and, con- 
versely, they should be taken into account in future experimental analyses aiming at extracting 
possible signals of new physics. 
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Introduction 



It has already been clearly established [l]-[28] that large Sudakov logarithms arising at TeV energy 
scales as a consequence of a non-cancellation between real and virtual contributions can enhance 
the effects of Electro- Weak (EW) corrections in electron-positron scattering, so that the latter grow 
as Og W log 2n (s/M^) at the n-th perturbative order even in fully inclusive observables, where s e + e - 
is the collider centre-of-mass (CM) energy squared and M\y the W boson mass. Eventually, they 
can even surpass the corrections generated in QCD: e.g., in the total hadronic cross section at 
-y/Sg+g- w 800 GeV and above. 

The reason for this is intimately related to the violation of the Bloch-Nordsieck theorem oc- 
curring in non-Abelian theories whenever the initial state has a finite (weak) isospin charge 5 , as 
dictated by the given beam configuration. This is immediately evident for leptonic colliders, as 
the Sudakov logarithms present in e + e~ scattering would cancel against those originating in e + u e 
and v e e~ collisions (the (anti)neutrinos are the isospin partners of the electron/positrons), a con- 
dition which is clearly impossible to satisfy experimentally. One can view the mechanism rather 
intuitively from a diagrammatic perspective. In short, virtual W corrections simply multiply the 
leading-order (LO) scattering matrix elements, thus being proportional to cr e + e - , while the real 
emission of a W boson does change the isospin of the incoming electron/positron and turns it into 
a(n) (anti)neutrino, so that the corrections here are proportional to (J e + Ve and 0"p ee -- 

Evidently, this does not occur for the case of real and virtual Z boson corrections (or photons, 
for that matters). The source of the large logarithms is then in principle manifest only in the 
case of W boson corrections. In practice, though, one should recall that both W and Z real 
bosons are unstable and decay into high transverse momentum leptons and/or jets, which are 
normally captured by the detectors. In the definition of an exclusive cross section, one may then 
remove events with such additional particles. Ultimately, other than being a second source of 
Bloch-Nordsieck violation for the case of W corrections, this merely experimental procedure will 
also spoil the cancellations between real and virtual contributions in the case of Z bosons, simply 
because the former are not included in the definition of the measured quantity. 

The leading, double-logarithmic, angular-independent weak logarithmic corrections are univer- 
sal, i.e. they depend only on the identities of the external particles. Both leading and subleading 
corrections are finite (unlike in QED and QCD), as the masses of the W and Z gauge bosons provide 
a physical cut-off for the otherwise divergent infrared behaviour. In some instances large cancel- 
lation between angular-independent and angular-dependent corrections [14] and between leading 
and subleading corrections [23] have been found at TeV energies. It is therefore of paramount 

5 The problem is in principle present also in QCD, with respect to the colour charge; in practice, however, it has 
no observable consequences, because of the final averaging of the colour degrees of freedom of the incoming partons, 
forced by their confinement into colourless hadrons. 
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importance to study the full set of fixed order weak corrections in order to establish the relative 
size of the different contributions at the energy scales which will be probed at TeV scale machines. 

Furthermore, weak contributions can be isolated in a gauge-invariant manner from purely 
Electro-Magnetic (EM) (or QED) effects [3], [7]-[9], which may or may not be included in the 
calculation, depending on the observable being studied and the aimed at accuracy. In view of all 
such arguments, it is then legitimate and topical to investigate the importance of higher-order weak 
effects at TeV scale hadronic colliders [14], such as Tevatron {\fs p p = 2 TeV) and LHC {\fs pp = 14 
TeV). 

Some further considerations are however in order in the hadronic context. First, one should 
recall that hadron-hadron scatterings (pp,pp) involve valence (or sea) partons of opposite isospin 
in the same process. Thus the above-mentioned cancellations may potentially be restored. For 
example, in pp(pp) scatterings one finds both uu(uu) and ud(ud) subprocess contributions to the 
total hadronic cross section, which tend to balance each other, this effect being actually modulated 
by the Parton Distribution Functions (PDFs). Secondly, several crossing symmetries among the 
involved partonic subprocesses can also easily lead to more cancellations. Thirdly, whether or not 
these two mechanisms take place, spin asymmetries due to weak effects would always be manifest 
in some observables, since QCD has a trivial helicity structure (just like QED). 

The purpose of this paper is that of establishing the importance of one-loop weak effects in 6-jet 
production at Tevatron and LHC. This is a pressing problem, as the pt distribution of Tevatron 
data for fo-quark production shows a clear disagreement with the theory [29], now known to next- 
to-leading order (NLO) accuracy in QCD [30] 6 , even after all uncertainties related to the definition 
of the cross section [31] and the extraction of the 6-quark fragmentation function are properly taken 
into account [32]. In order to avoid such uncertainties, we consider in this paper the cross section 
for di-jet production for which each jet contains a 6(6)-quark. Data from Run 2 is also expected to 
be presented in this format [29]. Comparisons of such 6-jet cross sections from Run 1 with NLO 
QCD [33] show a less severe discrepancy than in the case of 6-quark distributions. The comparison 
between theory and 6-quark/jet data is eventually expected to continue at LHC with much higher 
precision [34]. 

Production of 6-jets at Tevatron and LHC 

Even if the discrepancy referred to at the end of the previous section may not appear alarming at 
this stage, it is conceivable that the higher statistics available after Run 2 will afford the possibility 
of looking at more exclusive observables, in order to understand whether the difference may be due 
to some possible new physics effects, such as, e.g., W and Z' gauge bosons [35]. In this respect, 
it is natural to turn to quantities which are insensitive to QCD effects, such as the aforementioned 
6 Also the subleading LO tree-level contributions from EW interactions have been calculated. 
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spin induced asymmetries in the cross section. From this point of view, the knowledge of the weak 
effects described above would be of paramount importance, even if their overall contribution to the 
inclusive cross section should turn out to be negligible. 

After Run 2 at Tevatron, the accumulated statistics will be sufficient to select hadronic samples 
with two 6-jets and to establish their charge as well: e.g., by extracting two displaced vertices and 
measuring the charge of one of the (at least two) associated jets, via a high px lepton selection or 
jet charge reconstruction. This will enable one to define the usual 'forward-backward asymmetry' 
for 6-jets also at hadronic colliders, hereafter denoted by ^4fb 7 - Unfortunately, because of the 
symmetric beam configuration at LHC, one cannot define the forward-backward asymmetry in 
this case. Pure QCD contributions through orders a| and a| to such a quantity are negligible at 
Tevatron compared to the tree-level EW ones, which are of order a| w . We set out here to compute 
one- loop virtual effects to bb production through order a|aEW) which have then formally a similar 
strength to the purely EW ones, given that q| ~ oew at TeV energies. 

Before proceeding, we should like to clarify here that we will only include (in the language of 
Ref. [31]) 'flavour creation' contributions and neglect both the 'flavour excitation' and 'shower/frag- 
mentation' ones. While this is certainly not justified in the total inclusive 6-cross section [31], it 
is entirely appropriate for the bb one that we will be using in the definition of ^4fb, for which we 
will require 'two' high px 6-jets (thus depleting the 'flavour excitation' terms) tagged in opposite 
hemispheres (thus suppressing the 'shower/fragmentation' contributions). Finally, as anticipated in 
the previous discussion, we will neglect including QED corrections at this stage of our computation 
(this is indeed a gauge- invariant procedure, as we have explicitly verified), since we will ultimately be 
most interested in the forward-backward asymmetry, to which pure EM terms contribute negligibly. 

Partonic contributions to the pp/pp — > bb cross section 

The inclusive 6-jet cross section at both Tevatron and LHC is dominated by the pure QCD contri- 
butions gg — > bb and qq — > bb, known through order for n = 2, 3. Of particular relevance in this 
context is the fact that for the flavour creation mechanisms no as aw tree- level contributions are 
allowed, because of colour conservation: i.e., 



q b Q b 




where the wavy line represents a Z boson (or a photon) and the helical one a gluon. Tree-level 
asymmetric terms through the order a EW are however finite, as they are given by non-zero quark- 
7 In this respect, it is intriguing to recall the long-standing disagreement between data and SM for such a quantity, 
as seen at LEP and SLD [36], as well as several other observables involving &-quarks/jets, both at collider and fixed 
target experiments [37]. 
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antiquark initiated diagrams such as the one above wherein the gluon is replaced by a Z boson 
(or a photon). The latter are the leading contribution to the forward-backward asymmetry (more 
precisely, those graphs containing one or two Z bosons are, as those involving two photons are 
subleading in this case, even with respect to the pure QCD contributions). 

Here, we will compute one-loop and (gluon) radiative contributions through the order a^ay/, 
which - in the case of quark-antiquark induced subprocesses - are represented schematically by the 
following diagrams: 



X 

I 

X' 



b q i 




U + crossed box + 



1' + crossed box + 



]t + other three vertices + 



+ all self-energies + 



]' + gluon permutations. 

(2) 

The gluon bremsstrahlung graphs are needed in order to cancel the infinities arising in the virtual 
contributions when the intermediate gluon becomes infrared. Furthermore, one also has to include 
«| a w terms induced by gluon-gluon scattering, that is, interferences between the graphs displayed 
in Fig. 1 of Ref. [38] and the tree-level ones for gg — > bb. In the remainder of this paper, we will 
assume m& = and m t = 175 GeV (with Y t = 1.55 GeV): the top-quark enters the vertices and 
self-energies of the diagrams in (2) as well as the boxes (in additions to self-energies and vertices 
themselves) in Fig. 1 of Ref. [38], whenever a virtual W exchange occurs. The Z mass used was 
Mz = 91.19 GeV and was related to the W mass, M\v, via the SM formula Myy = Mzcosdw, 
where sin 2 6 W = 0.232. (Corresponding widths were T z = 2.5 GeV and T w = 2.08 GeV.) For a s 
we have used the one- or two-loop expressions as specified below, with A^- -4 -* set according to the 
PDFs used. 



q b 9 b 
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Some of the diagrams contain ultraviolet divergences. These have been subtracted using the 
'modified' Minimal Subtraction (MS) scheme at the scale fi = Mz- Thus the couplings are taken to 
be those relevant for such a subtraction: e.g., the EM coupling, «em = «ew sin 2 9\y, has been taken 
to be 1/128 at the above subtraction point. The one exception to this renormalisation scheme has 
been the case of the self-energy insertions on external fermion lines, which have been subtracted 
on mass-shell, so that the external fermion fields create or destroy particle states with the correct 
normalisation. 

Infrared divergences occur when the virtual or real (bremsstrahlung) gluon is either soft or 
collinear with the emitting parton. It is because we are considering 6-jets which include a possible 
gluon parallel to the 6-quark rather than open 6-quark production that the collinear divergences 
cancel, this way removing the logarithmic dependence on the 6-quark mass which was investigated 
and resummed in the analysis of Ref. [39]. Moreover, in our case the collinear divergences cancel 
amongst themselves. This can be seen since by colour conservation only interferences between 
gluon emission from the initial and final state quarks are permitted. If the gluon is parallel to 
an initial (final) quark then from the collinear vertex it is contracted into its own momentum and 
the sum of amplitudes for a longitudinal gluon emitted from both final (initial) states cancels by 
virtue of a Ward identity. For virtual corrections, the infrared divergences arise from the box 
graphs and there is an equivalent cancellation of collinear divergences between the crossed and 
uncrossed boxes. This leaves the soft divergences which can be readily extracted and as expected 
cancel between the virtual corrections and bremsstrahlung emissions. Nevertheless, for the sake 
of numerical stability when carrying out the necessary numerical integration over phase space and 
convolution with the PDFs, it is preferable to use the formalism of Catani and Seymour [40], 
whereby corresponding dipole terms are subtracted from the bremsstrahlung contributions in order 
to render the phase space integral free of infrared divergences. The integration over the gluon 
phase-space of these dipole terms are performed analytically in d— dimensions, yielding pole terms 
which cancel explicitly against the pole terms of the box graphs. 

Our expressions for each of the diagrams contain the complete helicity information from both 
the initial and final state. They have been calculated using FORM [41] and reproduced by an 
independent program based on FeynCalc [42]. The formulae have all been checked for gauge 
invariance. The full expressions for the contributions from all possible OgOEW graphs are too 
lengthy to be reproduced here. 

Numerical results for Tevatron and LHC 

We start our numerical investigation of the processes pp/pp — > bb by first computing the total cross 
section, a(pp — > bb), for Tevatron (Run 2). This can be found in Fig. 1 (top), as a function of the 
transverse momentum of the 6-jet (or 6-jet) and decomposed in terms of the various subprocesses 
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discussed so far. (Hereafter, the pseudorapidity is limited between —2 and 2 in the partonic CM 
frame.) The dominance at inclusive level of the pure QCD contributions is manifest, over the 
entire pr spectrum. At low transverse momentum it is the gluon-gluon induced subprocess that 
dominates, with the quark-antiquark one becoming the strongest one at large pr- The QCD K- 
factors, defined as the ratio of the a\ rates to the <x§ ones are rather large, of order 2 and positive 
for the gg — ► bb subprocess and somewhat smaller for the qq — ► bb case, which has a p^-dependent 
sign 8 . The tree-level a| w terms are much smaller than the QCD rates, typically by three orders 
of magnitude, with the exception of the px ~ Mz/2 region, where one can appreciate the onset 
of the Z resonance in s-channel. All above terms are positive. The a|aEW subprocesses display 
a more complicated structure, as their sign can change over the transverse momentum spectrum 
considered, and the behaviour is different in qq — > bb(g) from gg — ► bb. Overall, the rates for the 
a|aEW channels are smaller by a factor of four or so, compared to the tree-level a EW cross sections. 
Fig. 1 (bottom) shows the percentage contributions of the ag, a| w and a§aEW subprocesses, with 
respect to the leading «g ones, defined as the ratio of each of the former to the latter 9 . The a|aEW 
terms represent a correction of the order of the fraction of percent to the leading a| terms. Clearly, 
at inclusive level, the effects of the Sudakov logarithms are not large at Tevatron, this being mainly 
due to the fact that in the partonic scattering processes the hard scale involved is not much larger 
than the W and Z masses. 

Next, we study the above mentioned forward-backward asymmetry, defined as follows: 

^ = v+(PP -> bb) - a_(pp -> bb) 
a + {pp — ► bb) + a- (pp — ► bb) ' 

where the subscript +(— ) identifies events in which the 6-jet is produced with polar angle larger (smal- 
ler) than 90 degrees respect to one of the two beam directions (hereafter, we use the proton beam as 
positive z-axis). The polar angle is defined in the CM frame of the hard partonic scattering. Notice 
that we do not implement a jet algorithm, as we integrate over the entire phase space available 
to the gluon. In practice, this corresponds to summing over the two- and three-jet contributions 
that one would extract from the application of a jet definition. The solid curve in Fig. 2 (top) 
represents the sum of the tree-level contributions only, that is, those of order and a EW , whereas 
the dashed one also includes the higher-order ones a s and a§aEW- (Recall that the contributions 
to the asymmetry due to the pure QED and QCD terms a| M , a| and a| are negligible 10 .) 

The effects of the one-loop weak corrections on this observable are extremely large, as they are 
not only competitive with, if not larger than, the tree-level weak contributions, but also of opposite 

8 Further notice that in QCD at NLO one also has (anti)quark-gluon induced (tree-level) contributions, which are 
of similar strength to those via gluon-gluon and quark-antiquark scattering but which have not been shown here. 

9 In the case of the a| corrections, we have used the two-loop expression for as and a NLO fit for the PDFs, as 
opposed to the one-loop formula and LO set for the other processes (we adopted the GRV94 [43] PDFs with MS 
parameterisation) . 

10 And so would also be the one-loop a|aEM terms not computed here. 
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sign over most of the considered pt spectrum. In absolute terms, the asymmetry is of order —4% at 
the W, Z resonance and fractions of percent elsewhere, hence it should comfortably be measurable 
after the end of Run 2. 

Fig. 3 shows the same quantities as in Fig. 1, now defined at LHC energy. By a comparative 
reading, one may appreciate the following aspects. Firstly, the effects at LHC of the a§aEW 
corrections are much larger than the a EW ones already at inclusive level (see top of Fig. 3), as 
their absolute magnitude becomes of order —2% or so at large transverse momentum (see bottom 
of Fig. 3): clearly, logarithmic enhancements are at LHC much more effective than at Tevatron 
energy scales 11 . Secondly, the overall production rates at the CERN collider are in general much 
larger than those at FNAL, because of the much larger gluon component of the proton. 

Conclusions 

In summary, we should like to remark upon the following aspects of our analysis. 

• Inclusive corrections to the bb cross section due to one-loop weak interaction contributions 
through order a|aEW are small and undetectable at Tevatron, while becoming visible at 
LHC, because of the much larger cross section and luminosity available. In practice, the 
weak Sudakov logarithms are threshold suppressed at the FNAL collider while at the CERN 
machine they become sizable. In the former case then, they cannot explain the current data 
vs. theory discrepancy seen in the 6-quark/jet cross sections. 

• One-loop weak effects onto 6-quark asymmetries (e.g., we have studied the forward-backward 
one) are found to be large at Tevatron, where they can be defined experimentally. Here, the 
forward-backward asymmetry is subject to large corrections because the tree-level (quark- 
antiquark) subprocesses are formally of the same order as the one- loop contributions (initiated 
by both quark-antiquark and gluon-gluon collisions), eventually being measurable if collider 
luminosity plans will turn out to be on schedule. 

In conclusion, at both current and planned TeV scale hadronic colliders, one-loop weak effects from 
SM physics may be important and need to be taken into account particularly in order to extract 
possible signals of new physics from data. 
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Figure 1: The total cross section contributions vs. the transverse momentum of the 6-jet for pp — > bb 
production at Tevatron (2 TeV) as obtained via the various subprocesses discussed in the text (top) 
and the corrections due to the a^ w , a|a| w and a| terms relative to the a| ones (bottom). 
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Figure 2: The forward-backward asymmetry vs. the transverse momentum of the fr-jet for pp — > bb 
events at Tevatron (2 TeV), as obtained at tree- level and one-loop order (top) and the relative 
correction of the latter to the former (bottom). (Errors in the ratio are statistical.) 
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Figure 3: The total cross section contributions vs. the transverse momentum of the 6-jet for pp — > bb 
production at LHC (14 TeV) as obtained via the various subprocesses discussed in the text (top) 
and the corrections due to the a^ w and ag«|w terms relative to the a| ones (bottom). (Here, 
we do not show the corrections due to a| terms as results are perturbatively unreliable, given that 
iC-factors as large as 3-4 can appear.) 
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